Reentrant Behavior of the Spinodal Curve in a Nonequilibrium Ferromagnet 
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The metastable behavior of a kinetic Ising-like ferromagnetic model system in which a generic 
type of microscopic disorder induces nonequilibrium steady states is studied by computer simulation 
and a mean-field approach. We pay attention, in particular, to the spinodal curve or intrinsic 
coercive field that separates the metastable region from the unstable one. We find that, under 
strong nonequilibrium conditions, this exhibits reentrant behavior as a function of temperature. 
That is, metastability does not happen in this regime for both low and high temperatures, but 
instead emerges for intermediate temperature, as a consequence of the non-linear interplay between 
thermal and nonequilibrium fluctuations. We argue that this behavior, which is in contrast with 
equilibrium phenomenology and could occur in actual impure specimens, might be related to the 
presence of an effective multiplicative noise in the system. 
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I. INTRODUCTION 

The concept of metastability Q, Q is crucial to many 
branches of science. Metastable states occur in liq- 
uids and glasses, quark/gluon plasmas^ globular 
proteins, cosmological phase transitions, jjjthe "stan- 
dard model" of particle physics, Q climate models, Q 
black holes and protoneutronic stars, for instance. Un- 
derstanding metastability from a microscopic point of 
view is therefore most interesting. It is also a difficult 
task, given that the concept is a kinetic feature which 
is not described by the Gibbs ensemble theory. Con- 
sequently, a lot of activity still focuses on very simple 
cases, particularly, kinetic Ising-like models. Some recent 
studies along this line concern the details of nucleation 
during the relaxation processe s. 1 1(1 some exact results in 
the limit of zero temperature, [llj the checking of theo- 
retical predictions by means of computer simulation, [l^ 
and analysis of the effects of open borde rs. Il3j quenched 
impurities, an d demagnetizing fields. |l5| 

These studies deal with systems in which the 
metastable state decays towards the equilibrium sta- 
ble state. In this case, some understanding can be 
achieved via nucleation theories in which Gibbs ther- 
modynamic (equilibrium) free energy plays a central 
role. However, more general and intriguing is the case 
in which relaxation is towards a nonequilibrium steady 
state. 0, 0, 0, 0| Nonequilibrium conditions appear 
ubiquitously in nature, and they characterize the evo- 
lution of most real systems. [It) Under such conditions, 
no free energy can be defined in general, ^jj and no co- 
herent theoretical framework exists that accounts for the 



observed far-from-equilibrium behavior. Some important 
questions regarding metastability concern the existence 
and properties of a nonequilibrium macroscopic potential 
capturing the essential physics of the metastable-stable 
transition under nonequilibrium conditions, and the limit 
of metastability when such conditions hold. 

In this paper we therefore deal with aspects of metasta- 
bility in a kinetic Ising-like model with nonequilib- 
rium steady states. Our interest is on the effects of 
the nonequilibrium condition on the properties of the 
metastable state as one varies the system parameters. In 
particular, we study the magnetic-field strength for the 
onset of instability. This is the intrinsic coercive field |2(jj 
which locates in magnets the spinodal curve which is fa- 
miliar from studies of density-conserved systems. The 
system behavior around this curve is the consequence of 
a complex interplay between thermal and nonequilibrium 
fluctuations. This results in a spinodal curve that depicts 
novel behavior as compared to the equilibrium case. An 
interesting observation is that metastability occurs in the 
strong nonequilibrium regime at intermediate tempera- 
tures but not in the low temperature limit, pointing out 
that, in this regime, noise enhances metastability. 

Recent studies on critical behavior of some nonequi- 
librium models have predicted similar reentrance phe- 
nomena. That is, in a large class of model systems, one 
observes that, under nonequilibrium fluctuations, a dis- 
ordered phase which characterizes the system at both 
low and high temperatures becomes ordered at interme- 
diate temperatures. [H El |H |H E3 Further 
research is still needed, however, before one may con- 
clude on the relevance of such model behavior on the 
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reentrance phenomena reported, for instance, in nonequi- 
librium phase transitions driven by competition between 
quantum and thermal fluctuations in superconductors 
and vortex matter, [2^, Eij, l30| and concerning differ- 
ent li quid, g lassy and amorphous phases in water and 
silica, jj^J l32l l33| Despite the similarities between our re- 
sults and these studies, they are different in essence: the 
latter concern reentrance in phase diagrams associated to 
nonequilibrium phase transitions, while our work concern 
reentrance of the nonequilibrium spinodal curve, which 
characterizes the limit of metastability. 

The paper is organized as follows. Sections 11 and 111 
describe, respectively, the model and a dynamic mean- 
field approximation. Section IV contains our main results 
on the static properties of nonequilibrium metastability; 
in particular, we evaluate the intrinsic coercive field. In 
Section V we measure this spinodal field in computer 
simulations, and numerical results are compared in this 
section with our mean-field theory. Section VI is devoted 
to conclusions. 



II. THE MODEL 

Consider the two-dimensional Ising model on the 
square lattice of side L, A — {1, . . . , L} 2 C Z 2 , with peri- 
odic (toroidal) boundary conditions. There is a spin vari- 
able at each lattice site with two possible states, = ±1, 
isA, and spin-spin interactions and influence of an ex- 
ternal magnetic field h as described by the function 
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where s = {si} and the first sum runs over all nearest- 
neighbor pairs. Furthermore, the spin system evolves with 
time via stochastic single spin -flip dynamics as deter- 
mined by the master equation: 



dP(s;f) 
dt 



[w(s l -> s)P(sV) - w(s -> s l )P(s; t)] 



(2) 

Here, P(s; t) is the probability of configuration s at time 
t, s l stands for s after performing a flip at i, i.e., Sj — > — s,-, 
and lv(s — > s l ) stands for the corresponding transition 
rate. In order to ensure nonequilibrium conditions, we 
introduce a weighted competition between two differ- 
ent temperatures (one "infinite" and the other finite). 
This has been shown to be the simplest way of inducing 
nonequilibrium behavior in lattice models. |l7| The rate 
is then chosen to be 



uj(s -> s l ) = p + (I -p) 



i) 



1 + Q-p&Hisum) ' 



(3) 



where j3 — 1/T stands for the lattice (inverse) tem- 
perature — so that we are fixing the Boltzmann's con- 
stant to unity — , and AH(si,rii) = H(s l ) — H(s) = 
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FIG. 1: Two examples of spin domains, each following from a 
different lattice partition P(A); see the main text. The spins 
that do not belong to the domain are in black, surface spins 
are gray, and the spins forming the domain interior are empty. 



2si [2(rii — 2) + h], where n.j <E [0,4] is the number of up 
nearest-neighbors of s i . 

The rate © describes spin flips under the action of two 
competing heat baths. For p = 0, the system goes asymp- 
totically towards the unique Gibbs, equilibrium state for 
temperature T and energy 7i. This has a critical point 
at h = and T = T c (p = 0) = T ons « 2.2691, the On- 
sager, equilibrium critical temperature. For < p < I, 
the conflict in @ impedes canonical equilibrium, and 
the system evolves asymptotically towards a nonequi- 
librium steady state which may essentially differ from 
a Gibbs state. 0,13 I n this case, no equilibrium ther- 
modynamic global temperature can be defined. Now pa- 
rameter T can be thought as a source of thermal fluc- 
tuations, which compete with the non-thermal (nonequi- 
librium) noise induced by p. The system now exhibits a 
critical point, at h = and T = T c (p) < T r _ (0) , which 
is apparently of the Ising universality class, [33 |3g but 
only as far as p is small enough. In fact, the nonequi- 
librium disorder which is implied by washes out the 
critical point for any p > p c « O.17.0| O ne ma y think 
of the dynamic random perturbation parameterized by p 
as an extra source of (nonequilibrium) disorder and ran- 
domness which is likely to occur also in many actual sys- 
tems in nature. |16l Il7| A main question here is how the 
metastable states in the system depend on the compe- 
tition between this non-thermal noise and the standard 
thermal fluctuations parameterized by T. 



III. A MEAN-FIELD APPROXIMATION 

We first studya mean-field solution of in the pair 
approximation. |l7l l37| This approach is a dynamic gen- 
eralization of Kikuchi's method known as Cluster Varia- 
tion Method. [38| Consider a partition P of the lattice such 
that resulting domains, qj E P(A), satisfy qj n qy = 
if 3 7^ f, an d Uj Qj = A. We define the surface Sj of qj 
as the set of all its spins that have at least one nearest 
neighbor outside the domain; |3^| the rest is the interior, 
namely, Xj = qj — Sj , and it follows that qj = Xj U Sj . 
These definitions are illustrated in Fig^ Next, consider 
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a local observable A(s qj ; j) which exclusively depends on 
spins belonging to qj. One readily has from for the 
average (A{j)) t = Y^ s A(s qj ;j)P{s;t) that 



d(A(j)) 



-=EE &A(s qj ;j;i)u(s qj -> s\.)Q(s qj ;t) 



dt 

+ E E AA K ; *M S -> s *) p ( s ; *)> ( 4 ) 

s ieSj 

where s qj is the configuration of the domain spins, 
AA(s qj ;j;i) = A(s q .;j) - A(s qj ;j), and Q(s qj ;t) = 
J2 s -s P( s > ' s ^ ne probability of having the configu- 
ration s qj at time t. The notation uj(s qj — > s*.) in eq. 
Q stresses the fact that flipping a spin in the interior 
only depends on the spins belonging to the domain. 

Let us assume that the system is spatially homoge- 
neous, namely, that (A(j)) = (A), qj = q, Tj = T, and 
Sj = S for any j. Equivalently, the partition P(A) is 
regular, so that all domains are topologically identical. 
One notices that the two r.h.s. terms in eq.JU concern 
the domain interior and surface, respectively; the latter 
couples the domain dynamics to its surroundings. Our 
second approximation consists in neglecting this surface 
term, i.e., any correlation during time evolution which 
extends outside the domain. Under these two approxi- 
mations, homogeneity and kinetic isolation, equation 1|4|1 
reduces to 



d(A) t 
dt 



^^A4(s, ;i V( S , -> s q )Q(s q ;t). 



(5) 



Next, one needs to estimate Q(s q ;t) in terms of 
n— body correlation functions. Assuming that only (s) 
and (sjSj), with i and j nearest-neighbor sites inside 
the domain, matter, Q(s q ;t) may be written as a func- 
tion of the spin density p(s) and the density p(s, s') 
of nearest-neighbors pairs only. Furthermore, as only 
nearest-neighbors correlations are taken into account, we 
consider a domain with only one spin in its interior, which 
has 4 surface spins; see Fig^. With this choice, our 
mean-field theory turns out to be a nonequilibrium ana- 
log of the equilibrium Bcthe-Peierls approximation, ft 
follows that the probability of finding the central spin in 
state s surrounded by n up nearest-neighbor spins is 



Q(s q ;t) = Q{s,n) 



4W+, S )>(-, S ) 4 
n ) 



p{sf 



(6) 



Therefore, using the relations p(+, — ) = p(— ,+) = 
P{+) - P(+, +) and p{-, -) = 1 + p(+, +) - 2p(+), and 
writing x = p(+) and z = p(+, +), cq. J^J) reads 



d(A) t 



E 



dt £ — ' \ n 

71=0 



. (x-z) n (l + z~2x) 4 - n . 
-^A(-.n) ^_ g) 3 '—u,(-n) 



, (7) 
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FIG. 2: Variation with p of the critical temperature for the 
nonequilibrium ferromagnetic system in the first-order mean- 
field approximation. The solid line in the inset stands for the 
locally-stable steady magnetization as a function of temper- 
ature (in units of the equilibrium critical value) for h = 
and, from top to bottom, p = 0, 0.001, 0.005 and 0.01. The 
symbols in the inset are Monte Carlo results for a 53 x 53 
lattice. 



s*). This is for local isotropic 



where w{s,n) = lu(s q 

observables for which the dependence on s g is through 
the pair (s,n) only, A(s q ;t) = A{s,n). 

One may apply Q to the observables A\(s,ri) — 
\ (1 + s) and A2(s, n) = |n (1 + s) whose averages are x 
and z, respectively. Then, AA±(s,ri) — — s, &.Ai(s,tl) = 
— \sn, and 



d.r 

dt 

dz 
dt 



F 1 (x,z) = J2G(x,z;n) (8) 

4 

F 2 (x,z) = Yl jG(x,z;n), (9) 



71=0 



where 
G {x, z; n) 



(x - z) n (l + z- 2x) 4 



(1 - xf 
z n (x-z) 4 - n , 



-uj(-,n) 



Together with and (J7J), these equations provide x(t) 
and z(t) as well as any other local isotropic magnitude. 



IV. STATIC PROPERTIES 

Our main interest here is on the steady solutions 

F 1 (x st ,z st )=0, F 2 {x sU z st ) = 0, (10) 

and on their stability. [ioT] Both stable and metastable 
states are locally stable under small perturbations, which 
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requires the (necessary and sufficient) condition: |41j 



dF\ 

dx 



dFi 
dx 

m 

dz 



The condition 



dF 2 

dz 

dz 



dFi(x, z) 



< 0, 



and 



dF 2 
dx 



dx 



= 0, 



>0. (11) 



(12) 



on the other hand, characterizes incipient or marginal 
stability, i.e., the presence of a critical point {x° st , z^ t ) for 

|, and it follows that 



0, x 



si. 



\ and z c st 



T c (p) = -4 



1-Ap 



(13) 



This is the critical temperature for the nonequi- 
librium model in the present first-order mean-field 
approximation; [l7| see Fig[3 It is noticeable the exis- 
tence of a critical value of p such that T c (p c ) — 0, which 
gives p c — 5/32 = 0.15625 (to be compared with the 
exact value p c ~ 0.17). 

The stationary state (x s t, z s t) may be obtained numer- 
ically from the non-linear differential equations (|10fl sub- 
ject to the local stability condition (fill) - For h = 0, the 
up-down symmetry leads to pairs of locally-stable steady 
solutions, namely (x st ,z st ) and (1 - x st , 1 + z st - 2x st ). 
The result is illustrated in the inset of Fig. this 
also shows a comparison with Monte Carlo results which 
confirms the expected agreement at low and intermedi- 
ate temperatures for any p. The fact that increasing 
p at fixed T decreases the magnetization implies that 
the nonequilibrium perturbation tends to increase dis- 
order. For small enough fields, the situation closely re- 
sembles the case h = 0; the up-down symmetry is now 
broken, however, and locally-stable steady states with 
magnetization oriented opposite to the applied field are 
metastable. 

The locally-stable steady magnetization exhibits two 
branches as a function of h; see inset in Fig. [3J This 
hysteresis loop reveals that metastability does not occur 
for any \h\ > h* (T, v) > 0, where h*(T,p) is the intrinsic 
coercive Held. |42| Let z = z(x) the solution of eq. 
and write eq. (jHJ) as 



dec 
lit 



SV(x) 
Sx 



(14) 



This defines V(x), a (nonequilibrium) bimodal potential 
that controls the system time evolution. An increase in 
the field tends to attenuate the local minimum in V(x) 
associated with the metastable state. This minimum ex- 
ists only for \h\ < h*(T,p); the set of eas. HlOJI has only 
one solution, with magnetization sign equal to that of the 
applied field, for \h\ > h*(T,p). 




Ill 



FIG. 



3: h*(T,p), as a function of T for, from top to bottom, 
p = 0, 0.01, 0.02, 0.03, 0.031, 0.032, 0.035, 0.04, 0.05 and 
0.1. The qualitative change of behavior in the low temper- 
ature region occurs for p 6 (0.031,0.032). Inset: The two 
locally-stable steady magnetization branches as a function of 
h for T = 0.7T C (0) and p = 0.005. The solid (dashed) line 
represents stable (metastable) states. The dot-dashed line 
signals the discontinuous transition, at h — h*(T,p), where 
metastable states disappear. 



In order to evaluate h*(T,p), one may study how the 
metastable state responds to small perturbations of the 
applied field. If ^^z^j is a locally-stable stationary 

state for T, p and ho, with magnetization opposed to 
h, and we perturb h = ho + Sh, the new locally-stable 



stationary solution is modified: x st 



~h 



e, . One obtains at first order that 



dF 2 dF\ 
dh dz 



dF\ dF 2 
dh dz 



dx dz 



dF 2 dFi 
dx dz 



and z h st 



Sh, (15) 



flQ HQ , _ 



This (linear) response diverges for 



dFi dF 2 dF 2 dF 1 



dx dz 



dx dz 



X^°,Z^,h ;T, P 



0, (16) 



which corresponds to a discontinuity in the metastable 
magnetization as a function of h. For fixed T and p, 
the field for which lO holds is h*(T,p). Fig. shows 
the mean-field result for h*(T,p). In particular, for 
p = (the equilibrium case) we recover the standard 
equilibrium mean-field spinodal curve: converging to 2 
as T — > 0, linearly decreasing with temperature for small 
T, and vanishing as (T c — T) 3 / 2 at the mean-field equi- 
librium critical point. A main result derived from Fig. 
[3] is the existence of two different low temperature lim- 
its for h*(T,p). For small enough values of p, namely, 
p 6 [0, 0.031], which includes the equilibrium case, p = 0, 
the field h*(T,p) monotonously grows and extrapolates 
to 2 as T — > 0. For larger p, namely, p G [0.032, ^), how- 
ever, h*(T,p) — > as T — > 0, exhibiting a maximum at 
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Class 


Central spin 


Number of up neighbors 


AH 


1 


+ 1 


4 


8J+2h 


2 


+1 


3 


4J+2h 


3 


+1 


2 


2h 


4 


+1 


1 


-4J+2h 


5 


+1 





-8J+2h 


6 


-1 


4 


-8J-2h 


7 


-1 


3 


-4J-2h 


8 


-1 


2 


-2h 


9 


-1 


1 


4J-2h 


10 


-1 





8J-2h 



TABLE I: Spin classes for the two-dimensional Ising model 
with periodic boundary conditions. The last column shows 
the energy cost of flipping the central spin. 



intermediate T. The value p = tt c w 0.0315 separates the 
two types of asymptotic behavior. 

When we cool the system in the regime p < ir c , the field 
h*(T,p) increases, so in this case we require a stronger 
field to destroy the metastable state. This may be un- 
derstood on simple grounds. The tendency of spins to 
line up in the direction of the field competes with the 
tendency to maintain order implied by their mutual in- 
teractions. A metastable state lasts for a long time be- 
cause the latter prevails over the action of the field. Both 
T and p induce disorder; therefore, lowering T increases 
order, so that a stronger field is needed to destroy the 
metastable state, which is in fact observed for p < ir c . On 
the other hand, as p is increased, the disorder increases, 
and h* (T, p) needs to decrease for a fixed T, according to 
our observations. 

The picture for p > n c is more intriguing. Consider 
the case \h\ — 0.25 and p — 0.05 > ir c . As illustrated in 
FigEl one may define two temperatures, T\ < T2, such 
that metastable states only occur for T e (Ti,T 2 ). The 
fact that h*(T,p) extrapolates to zero in the low temper- 
ature limit for p = 0.05 > ir c indicates that such amount 
of nonequilibrium noise is able to destroy on its own any 
metastability. Following the above reasoning, increasing 
T adds disorder, so that no metastability should, in prin- 
ciple, show up in this case. However, there is a regime 
of intermediate temperatures, T £ (Ti,T2), for which 
metastability occurs. This noise-enhanced metastability 
is a consequence of the complex interplay between the 
standard thermal fluctuations and nonequilibrium noise: 
although both noises add independently disorder, their 
combination determines the existence of regions in the 
parameter space (T,p) in which no metastable states 
occur at low T but only at intermediate temperatures. 
This reentrance phenomenon is reminiscent of the one ob- 
served in the annealed I sing m odel l22l and other closely 
related systems [H |H HlHl |H [Zfl where multiplica- 
tive noise seems to play an essential role. [2^ 



V. MONTE CARLO SIMULATIONS: GROWTH 
AND SHRINKAGE OF THE STABLE PHASE 

In this section, we check further our theoretical predic- 
tions against computer simulation data. With this aim, 
we need a simple criterion to conclude that the model 
system exhibits metastable states. Let us first charac- 
terize all the possible local configurations in terms of the 
spin state, s = ±1, and the number of its up nearest 
neighbors, n G [0,4]. For periodic boundary conditions, 
there are 10 different spin classes, as shown in Table [I] 
The cost AH(s, n) of flipping any spin within a class is 
the same. That is, the rate Q only depends on s and n, 
which define the class. If rife(s) is the number of spins in 
class k when the system is is configuration s, and notic- 
ing that classes k 6 [1,5] are characterized by a central 
up spin, we may write the number of up spins which flip 
per unit time in the state s as 



G(s) = ^n fc (s)w fc . 



(17) 



fc=i 



As far as h < 0, this is the growth rate of the stable phase 
in state s. The shrinkage rate of the stable phase follows 
similarly as, 01 



S(s) 



in 

E 

k=6 



n k (s)uj k . 



(18) 



Given a phase-space point s, the rates G(s) and S(s) 
yield the local tendency of the system to evolve toward 
the stable or metastable phases, respectively. 

For h < 0, a state with all spins up, Si = {s, = +1, i = 
l,...,N = L 2 }, will relax after some time toward the 
stable steady state, which corresponds in this case to a 
configuration with negative magnetization, m < 0. For 
a given experiment j of a total of N exp runs, this relax- 
ation will proceed through certain path in phase space, 



(j) 

2 ! ' 



Sp'^}. Here s 



(j) 



is 



which we note as Uj = {si,s 
the Z-th configuration, starting from Si, of a total number 
of r(j) configurations which make up the path in exper- 
iment j. At any stage s/^ of this path, the difference 
G(sP) — S'(sp' 1 ) defines the net tendency of the system to 
evolve toward the final steady stable state. A metastable 
state is characterized by the presence of free energy bar- 
riers hampering the relaxation toward the truly stable 
state. In this case, relaxation is an activated process 
controled by large fluctuations. On the other hand, an 
unstable state decays without any hindrance. Therefore, 
we may divide relaxation paths in two different types. 
On one hand, metastable-like paths, in which at least one 
configuration s£p S <7j exists, excluding the final stable 
one, such that G(s£p) — S'(s^' 1 ) < 0, and on the other 



hand, unstable-like paths, such that G(s^) 



s(s£ ] ) > 



Vs^' e dj, excluding again the final stable state. 



6 




FIG. 4: Monte Carlo results for h* (T, p) as a function of T 
for L = 53 and, from top to bottom, p = 0, 0.01, 0.03, 0.0305, 
0.0320, 0.0350, 0.04 and 0.05. Notice the change of asymptotic 
behavior in the low temperature limit for p G (0.03,0.0305). 
Inset: The probability of the metastable state, as defined in 
the main text, as a function of h < for L — 53, T = 0.7T C (0) 
and p = 0. Data here correspond to an average over 500 in- 
dependent demagnetization experiments for each value of h. 
Error bars are smaller than the symbol sizes. 



magnetization is to. One may then define 
G(m) = ^( S ) G ( S )> S(m) = V ( S ) S ( S )> ( 19 ) 

{s|m} {s|m} 

where {s|to} are all system configurations with fixed 
magnetization to, and V(s) is the probability of observ- 
ing a configuration s during the relaxation from the initial 
state, si, toward the stable one. Steady states are usu- 
ally determined by the condition G{m) = S(m). How- 
ever, the lack of intersection between the curves G(m) 
and S(m) for h < in the to > interval does not 
contain information about the limit between metastable 
and unstable states. |^ Instead, the magnetic field for 
which such lack of intersection first develops defines the 
so-called Dynamic Spinodal field, }idsp(T,p), which di- 
vides the metastable region of parameter space (T, p, h) 
for finite systems in two subregions characterized by 
different relaxation morphologies. |45| In particular, for 
\h\ < hosp{T,p) the metastable state relax through the 
nucleation of a single droplet of the stable phase, while 
for hosp{T,p) < \h\ < h*(T 7 p) the relaxation proceeds 
via the nucleation of multiple stable-phase droplets. 



For fixed T, p and h < 0, given the stochastic- 
ity of the dynamics, one needs to be concerned with 
the probability of occurrence of metastability, defined as 
n me t{T,p,h)/N exp , where n met (T,p,h) is the number of 
experiments out of the total N exp in which the relaxation 
path in phase space belongs to the class of metastable-like 
paths. This is shown in the inset of Fig@| The intrinsic 
coercive field, h*(T,p), is defined in this scheme as the 
field for which n met (T,p,h*)/N exp = 0.5; this is shown 
in Fig. 2] for a system with L = 53. 

A detailed comparison of these numerical results with 
the theory in Fig[3]depicts semi-quantitative agreement, 
namely, the agreement is excellent except — as one should 
have expected — near the critical temperature. In par- 
ticular, the numerical critical value ir c for p is tt^ ic ~ 
0.03025, rather close to the theoretical prediction 7r£ air w 
0.0315. This nicely confirm that the addition of sufficient 
thermal noise in the presence of a large enough nonequi- 
librium perturbation, p > 7r c , tends to restore metasta- 
bility. We have also looked for finite-size corrections to 
the measured spinodal field by simulating larger systems, 
finding that these corrections are very small, and can be 
neglected for all practical purposes. 

Finally, let us remark that the stable phase growth 
and shrinkage rates have been introduced in literature as 
projected on the slow observable characterizing the relax- 
ation process, namely the system magnetization to. 43] In 
this case, the rate G{m) (5(m)) yields the average num- 
ber of up (down) spins which flip per unit time when 



VI. CONCLUSION 

This paper deals with some of the static properties 
of metastable states in a nonequilibrium Ising-like fer- 
romagnetic model system, as obtained from first-order 
mean-field theory and computer (Monte Carlo) simula- 
tions. We studied, in particular, the spinodal or intrin- 
sic coercive field, h*(T,p), defined as the magnetic field 
strength for which the metastable state becomes unsta- 
ble. Our theoretical approximation predicts reentrace 
phenomena as a function of T in the strong nonequi- 
librium regime, p > n c w 0.0315, where p controls the 
dynamic, non-equilibrium perturbation. More specifi- 
cally, metastability is not observed at low temperatures 
for p > 7r c , but it occurs as one increases the temperature. 
This noise-enhanced metastability reveals the existence 
of a complex interplay between the thermal and nonequi- 
librium noises. That is, adding the two effects — which, 
independently, tend to increase disorder — not always re- 
sults in decreasing the system ordering. The above is 
fully confirmed in computer simulations, in which h* may 
accurately be measured from the stable phase growth and 
shrinkage rates. 

The physical origin of the observed noise-driven 
metastability is intriguing. A clue to understand this 
behavior is to notice that certain systems under the 
action of a multiplicative noise exhibit a similar reen- 
trant behavior, namely, disorder is dominant at the low 
and high temperature regimes while well-defined order 
sets in at intermediate temperatures. That the com- 
petition between thermal and nonequilibrium fluctua- 
tions in may induce an effective multiplicative noise 
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can be understood on simple grounds. The effect of 
the nonequilibrium perturbation in our model may be 
described by means of an effective temperature T e ff, 
which is inhomogeneous throughout the system for any 
p > 0. 19J In fact, for any AH ^ we may write © 
as an equilibrium Glauber rate with effective parameters, 
w = exp(-/3 e// A7i)/[l + exp(-0 rff AH)\, and this de- 
fines an effective temperature |4y] 



T e ff{s,n) 



AH(s,n) 



In [ 



(20) 



As a matter of fact, the temperature a spin effectively 
feels depends on the local order around it, i.e., on the 
number of nearest neighbors pointing in the same direc- 
tion; T e jf is an increasing function of local order and, 
consequently, the amplitude of the fluctuations depends 
on the local order parameter. This is a main feature 
of Langevin-type models with a multiplicative noise. [25) 
Developing further this possibility to treat the limit of 
metastability is an open question. This work, which 
seems most interesting, is outside the scope of present 
paper. 
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